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Abstract 

The paper discusses and surveys some aspects of the potential theory of subordinate Brown- 
ian motion under the assumption that the Laplace exponent of the corresponding subordinator 
is comparable to a regularly varying function at infinity. This extends some results previously 
obtained under stronger conditions. 
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1 Introduction 

An Revalued process X = {Xt : t > 0) is called a Levy process in W 1 if it is a right continuous 
process with left limits and if, for every s, t > 0, X t+S — X s is independent of {X r ,r 6 [0, s]} and 
has the same distribution as X s — Xq. A Levy process is completely characterized by its Levy 
exponent 3> via 

E[exp{i<£, X t - X )}] = exp{-t$(0}, i > 0,£ G M d 
The Levy exponent $ of a Levy process is given by the Levy-Khintchine formula 

HO = + M T ) + J (l- + i(£, z)l {W<1} ) U(dx), ££R d , 

where I G M. d , A is a nonnegative definite d x d matrix, and II is a measure on W 1 \ {0} satisfying 
/(l A |x| 2 )II(dx) < oo. A is called the diffusion matrix, IT the Levy measure, and (l,A,TL) the 
generating triplet of the process. 

Nowadays Levy processes are widely used in various fields, such as mathematical finance, actu- 
arial mathematics and mathematical physics. However, general Levy processes are not very easy 
to deal with. 

A subordinate Brownian motion in R d is a Levy process which can be obtained by replacing 
the time of Brownian motion in R d by an independent subordinator (i.e., an increasing Levy 
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process starting from 0). More precisely, let B = (Bt : t > 0) be a Brownian motion in M. d and 
S = (St : t > 0) be a subordinator independent of B. The process X = (X t : t > 0) denned 
by X t = i?5 t is a rotationally invariant Levy process in Mr and is called a subordinate Brownian 
motion. 

The subordinator S used to define the subordinate Brownian motion X can be interpreted as 
"operational" time or "intrinsic" time. For this reason, subordinate Brownian motions have been 
used in mathematical finance and other applied fields. Subordinate Brownian motions form a very 
large class of Levy processes. Nonetheless, compared with general Levy processes, subordinate 
Brownian motions are much more tractable. If we take the Brownian motion B as given, then 
X is completely determined by the subordinator S. Hence, one can deduce properties of X from 
properties of the subordinator S. On the analytic level this translates to the following: Let <f> denote 
the Laplace exponent of the subordinator S, that is, E[exp{— XSt}] = exp{— tcfi(X)}, A > 0. Then 
the characteristic exponent of the subordinate Brownian motion X takes on the very simple form 
= 0(|a;| 2 ) (our Brownian motion B runs at twice the usual speed). Therefore, properties of 
X should follow from properties of the Laplace exponent (f). 

The Laplace exponent ^ of a subordinator S is a Bernstein function, hence it has a representation 
of the form 



where b > and \i is a measure on (0, oo) satisfying ^(1 A t) fj,(dt) < oo. If fi has a completely 
monotone density, the function (j) is called a complete Bernstein function. The purpose of this work 
is to study the potential theory of subordinate Brownian motion under the assumption that the 
Laplace exponent eft of the subordinator is a complete Bernstein function comparable to a regularly 
varying functions at infinity. More precisely, we will assume that there exist a E (0, 2) and a 
function £ slowly varying at infinity such that 



Here and later, for two functions / and g we write /(A) X g(X) as A — > oo if the quotient f(X)/g(X) 
stays bounded between two positive constants as A — > oo. 

A lot of progress has been made in recent years in the study of the potential theory of subordinate 
Brownian motions, see, for instance p3J EJ EU [26j ETJ EU [35] and (6J Chapter 5]. In particular, an 
extensive survey of results obtained before 2007 is given in [6l Chapter 5]. At that time, the focus 
was on the potential theory of the process X in the whole M d , the results for (killed) subordinate 
Brownian motions in an open subset still being out of reach. In the last few years significant 
progress has been made in studying the potential theory of subordinate Brownian motions killed 
upon exiting an open subset of Mr. The main results include the boundary Harnack principle 
and sharp Green function estimates. For processes having a continuous component see [26J (for 
the one-dimensional case) and [121 [T3| HH] (for multi-dimensional case) . For purely discontinuous 
processes, the boundary Harnack principle was obtained in [25] and sharp Green function estimates 




(0,oo) 




0(A) x A Q / 2 £(A) , A -»• oo . 
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were discussed in the recent preprint [27J. The main assumption in [T3j[TU[25] an d [3 Chapter 5] 
is that the Laplace exponent of the subordinator is regularly varying at infinity. The results were 
established under different assumptions, some of which turned out to be too strong and some even 
redundant. Time is now ripe to put some of the recent progress under one unified setup and to 
give a survey of some of these results. The survey builds upon the work done in (6j Chapter 5] and 
|25j . The setup we are going to assume is more general than all these of the previous papers, so in 
this sense, most of the results contained in this paper are extensions of the existing ones. 

In Section [2] we first recall some basic facts about subordinators, Bernstein functions and com- 
plete Bernstein functions. Then we establish asymptotic behaviors, near the origin, of the potential 
density and Levy density of subordinators. 

In Section [3] we establish the asymptotic behaviors, near the origin, of the Green function and 
the Levy density of our subordinate Brownian motion. These results follow from the asymptotic 
behaviors, near the origin, of the potential density and Levy density of the subordinator. 

In Section H] we prove that the Harnack inequality and the boundary Harnack principle hold for 
our subordinate Brownian motions. 

The materials covered in this paper by no means include all that can be said about the potential 
theory of subordinate Brownian motions. One of the omissions is the sharp Green function estimates 
of (killed) subordinate Brownian motions in bounded C 1 ' open sets obtained in the recent preprint 
[27] . The present paper builds up the framework for [27] and can be regarded as a preparation for 
[27] in this sense. Another omission is the Dirichlet heat kernel estimates of subordinate Brownian 
motions in smooth open sets recently established in [HI [9], [THl [11] . One of the reasons we do not 
include these recent results in this paper is that all these heat kernel estimates are for particular 
subordinate Brownian motions only and are not yet established in the general case. A third notable 
omission is the spectral theory for killed subordinate Brownian motions developed in [17 1 118 ^ IT9]. 
Some of these results have been summarized in [341 Section 12.3]. A fourth notable omission is 
the potential theory of subordinate killed Brownian motions developed in [211 1121 EH ESI HOI HI]- 
Some of these results have been summarized in [6j Section 5.5] and [341 Chapter 13]. In this paper 
we concentrate on subordinate Brownian motions without diffusion components and therefore this 
paper does not include results from[l3l [HI [26]. One of the reasons for this is that subordinate 
Brownian motions with diffusion components require a different treatment. 

We end this introduction with few words on the notations. For functions / and g we write 
f(t) ~ g(t) as t — > 0+ (resp. t — > oo) if the quotient f(t)/g(t) converges to 1 as t — > 0+ (resp. 
t — > oo), and f(t) x g(t) as t — > 0+ (resp. t — > oo) if the quotient f(t)/g(t) stays bounded between 
two positive constants as t — > 0+ (resp. t — > oo). 
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2 Subordinators 



2.1 Subordinators and Bernstein functions 

Let S = (St : t > 0) be a subordinate-! - , that is, an increasing Levy process taking values in [0, oo) 
with So = 0. A subordinator S is completely characterized by its Laplace exponent 4> via 

E[exp(-ASt)] = exp(-i0(A)) , A > 0. 

The Laplace exponent <ft can be written in the form (cf. [21 p. 72]) 

[•OO 

0(A) = bX + (1 - e~ xt ) fi(dt) . 
Jo 

Here b > 0, and is a <r- finite measure on (0, oo) satisfying 

/•oo 

/ (t A 1) /ti(di) < oo . 

The constant 6 is called the drift, and \i the Levy measure of the subordinator S. 

A C°° function </> : (0, oo) — > [0, oo) is called a Bernstein function if (— l) n D n (p < for every 
positive integer n. Every Bernstein function has a representation (cf. |3lj Theorem 3.2]) 

0(A) = a + b\ + [ (1 - e- xt ) fi(dt) 

J(0,oo) 

where a, b > and \i is a measure on (0, oo) satisfying Jj Q ^(1 At) \i{d£) < oo. a is called the killing 
coefficient, b the drift and fi the Levy measure of the Bernstein function. Thus a nonnegative 
function 4> ° n (0, oo) is the Laplace exponent of a subordinator if and only if it is a Bernstein 
function with </>(0+) = 0. 

Sometimes we need to deal with killed subordinators, that is, subordinators killed at independent 
exponential times. Let e a be an exponential random variable with parameter a > 0, i.e., P(e a > 
t) = e~~ at , t > 0. We allow a = in which case e a = oo. Assume that S is a subordinator with 
Laplace exponent 4> and e a is independent of 5. We define a process S by 

s - I 5 *' * < 6a 

[oo t > e a 

The process S is the subordinator 5 killed at an independent exponential time. We call S a killed 
subordinator. The corresponding Laplace exponent cj) is related to 4> as 

0(A) = a + 0(A), A>0. 

In fact, 

E[e^ §t ] = E[e^ 5t l {t<ea} ] = E[e^ 5 ']P(i < e a ) = e ~^ W e - at = e -*( a +<^W). 
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A function : (0, oo) —> [0, oo) is the Laplace exponent of a killed subordinator if and only if is 
a Bernstein function. For this reason, we use a killed subordinator sometimes. 

A Bernstein function is called a complete Bernstein function if the Levy measure /i has a 
completely monotone density n(t), i.e., (—\) n D n [i > for every non-negative integer n. Here and 
below, by abuse of notation we will denote the Levy density by n(t). Complete Bernstein functions 
form a large subclass of Bernstein functions. Most of the familiar Bernstein functions are complete 
Bernstein functions. See [341 Chapter 15] for an extensive table of complete Bernstein functions. 
Here are some examples of complete Bernstein functions: 



(i) 0(A) = A"/ 2 , a e (0,2]; 

(ii) 0(A) = (A + m 2 l a ) a l 2 - m, a £ (0, 2), m > 0; 

(iii) 0(A) = \ a ' 2 + A^ 2 , < p < a € (0, 2]; 

(iv) 0(A) = A a / 2 (log(l + A))T/ 2 , a e (0, 2), 7 e (0, 2 - a); 

(v) 0(A) = A Q / 2 (log(l + A))-' 3 / 2 , < /3 < a € (0, 2]. 



An example of a Bernstein function which is not a complete Bernstein function is 1 — e . 

It is known (cf. [341 Proposition 7.1]) that is a complete Bernstein function if and only if 
the function A/0(A) is a complete Bernstein function. For other properties of complete Bernstein 
functions we refer the readers to [34J. 

The following result, which will play an important role later, says that the Levy density of a 
complete Bernstein function cannot decrease too fast in the following sense. 

Lemma 2.1 ([27, Lemma 2.1]) Suppose that is a complete Bernstein function with Levy den- 
sity [i. Then there exists C\ > such that ji(t) < Ciji(t + 1) for every t > 1. 

Proof. Since [i is a completely monotone function, by Bernstein's theorem (cf. |34} Theorem 1.4]), 
there exists a measure m on [0, oo) such that fj,(t) = Jj ^ e~ tx m(dx). Choose r > such that 
l\0r] e ~ X m (.dx) > Jj roo ) e~ x m(dx). Then, for any t > 1, we have 




> 



e 



(t-l)r 



(r,oo) 



e x m(dx) > 



L 



(r,oo) 



e m(dx). 



Therefore, for any t > 1 



A*(* + 1) > 



/ 

.7 re 



[0,r] 



e 



(t+i> 



m{dx) > e 



[0,r] 



e te m[dx) 
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□ 

The potential measure of the (possibly killed) subordinator S is defined by 

POD 

U(A)=E l {St eA}dt, Ac[0,oo). (2.1) 
J o 

Note that U(A) is the expected time the subordinator S spends in the set A. The Laplace transform 
of the measure U is given by 

roo roo -i 

CU{\) = / e~ xt dUit) = E / exp(-AS't) dt = — — . (2.2) 
Jo Jo <A(A) 

We call a subordinator S a complete subordinator if its Laplace exponent is a complete 
Bernstein function. The following characterization of complete subordinators is due to [411 Remark 
2.2] (see also [6, Corollary 5.3]). 

Proposition 2.2 Let S be a subordinator with Laplace exponent (f> and potential measure U. Then 
(f> is a complete Bernstein function if and only if 

U{dt) = c5 (dt) +u(t)dt 

for some c > and completely monotone function u. 

In case the constant c in the proposition above is equal to zero, we will call u the potential 
density of the subordinator S. 

An inspection of the argument, given in [6l Chapter 5] or [41] : leading to the proposition above 
yields the following two results (cf. jfH Corollary 5.4 and Corollary 5.5] or [41 \ Corollary 2.3 and 
Corollary 2.4]). 

Corollary 2.3 Suppose that S = (St : t > 0) is a subordinator whose Laplace exponent 



oo 



0(A) = b\+ (1 - e~ M ) n(dt) 
Jo 

is a complete Bernstein function with b > or /i(0, oo) = oo. Then the potential measure U of S 
has a completely monotone density u. 

Proof. By [6l Corollary 5.4] or [41 [ Corollary 2.3], if the drift of the complete subordinator S is 
zero or the Levy measure \x has infinite mass, then the constant c in Proposition 12.21 is equal to zero 
so the potential measure U of S has a density u. The completeness of the density follows directly 
from Proposition 12.21 □ 
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Corollary 2.4 Let S be a complete subordinator with Laplace exponent 0(A) = J^°(l — e~ xt )^(dt). 
Suppose that the Levy measure fi has infinite mass. Then the potential measure of a (killed) subor- 
dinator with Laplace exponent ip(X) := A/0(A) has a completely monotone density v given by 

v(t) = fi(t, oo). 

Proof. Since the drift of S is zero and the Levy measure /i has infinite mass, by [6, Corollary 5.5] 
or [HI Corollary 2.4], we have that 

poo 

■0(A) = a+ (l-e- xt )v{dt) 
Jo 

where a = (f °° tjji{t)dt) , the Levy measure v of tp has infinite mass and the potential measure 
of a possibly killed (i.e., a > 0) subordinator with Laplace exponent ifi has a density v given by 
v(t) = fj,(t, oo). The completeness of the density follows from [BJ Corollary 5.3], which works for 
killed subordinators. □ 



2.2 Asymptotic behavior of the potential and Levy densities 

From now on we will always assume that S is a complete subordinator without drift and that the 
Laplace exponent of S satisfies lim^oo 0(A) = oo (or equivalently, the Levy measure of S has 
infinite mass). Under this assumption, the potential measure U of S has a completely monotone 
density u (cf. Corollary 12 . 3[) . The main purpose of this subsection is to determine the asymptotic 
behaviors of u and (i near the origin. For this purpose, we will need the following result due to 
Zahle (cf. [SI Theorem 7]). 

Proposition 2.5 Suppose that w is a completely monotone function given by 

roo 

w(t)= / e- st f(s)ds, 
Jo 

where f is a nonnegative decreasing function. Then 

/(«) < (1 -e-y^^fs' 1 ), s>0. 
Lf, furthermore, there exist 5 6 (0, 1) and a, sq > such that 

w{\t) < a\~ 5 w(t), A>l,t>l/s , (2-3) 
then there exists C2 = C^iu, /, a, sq, S) > such that 

f(s) > C 2 s- 1 w(s- 1 ), s < s . 
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Proof. Using the assumption that / is a nonnegative decreasing function, we get that, for any 
r > 0, we have 



> - 



i r _ t ./s\ , . i 



e- s f( 7 )d S > ±/ M(l-e-). 



i ./n V t ) ~ V V t 

Thus 

„ /r\ twit) 

In particular, we have 

/(a)< (l-e- 1 )" 1 ^ 1 ^" 1 ), s >0. 

and 

/ (£) ^ ( X " ^T 1 ' s>0,t>0. (2.4) 

On the other hand, for r € (0, 1], we have 



= / e -'f[-)d8+ J e- S f{- t )ds 
1 I r-V(-)^ + /'''' 



where in the last line we used (12. 4|) . Now we assume (|2.3p . then we get that 

w {^~ \ — as 5 w(t), t > l/so,s < r. 

Thus, for r S (0, 1], we have, 

tw(t) < a (1 - e- 1 )' 1 tw(t) J e- s s 5 - 1 ds + f (£) e _r . 

Choosing r £ (0, 1] small enough so that 

a(l-e- 1 )" 1 f e^s^dsK -, 
Jo 2 

we conclude that for this choice of r, we have 

/ Q > C!tw(t), t > 1/SQ 

for some constant c\ > 0. Since w is decreasing, we have 

/(s) > cx-io ( - ) > c 2 s' 1 w(s~ 1 ), s<rs , 
s \s/ 



where C2 = c\r. From this we immediately get that there exists C3 > such that 

f(s) > C3S~ 1 w(s~ 1 ), s < s - 

□ 

Corollary 2.6 The potential density u of S satisfies 

u{t) < Cst- 1 ^- 1 )- 1 , t > . (2.5) 
Proof. Apply the first part of Proposition 12.51 to the function 

w{t) := / e~ st u(s) ds 







□ 



We introduce now the main assumption on our Laplace exponent <f> of the complete subordinator 
S that we will use throughout the rest of the paper. Recall that a function I : (0, 00) — > (0, 00) is 
slowly varying at infinity if 

lim [, ) = 1 , for every A > . 

t^oc e(t) 

Assumption (H): There exist a £ (0,2) and a function £ : (0, 00) — > (0, 00) which is measurable, 
locally bounded above and below by positive constants, and slowly varying at infinity such that 

</>(X) x \ a/2 £{\) , A -> 00 . (2.6) 



Remark 2.7 The precise interpretation of (|2.6|) will be as follows: There exists a positive constant 
c > 1 such that 

The choice of the interval [l,oo) is, of course, arbitrary. Any interval [a, 00) would do, but with 
a different constant. This follows from the continuity of <p an d the assumption that I is locally 
bounded above and below by positive constants. Moreover, by choosing a > large enough, we 
could dispense with the local boundedness assumption. Indeed, by [U Lemma 1.3.2], every slowly 
varying function at infinity is locally bounded on [a, 00) for a large enough. 

Although the choice of interval [l,oo) is arbitrary, it will have as a consequence the fact that 
all relations of the type f(t) x g{t) as t — > 00 (respectively t — > 0+) following from (12. 6|) will be 
interpreted as cT l < f(t)/g(t) < c for t > 1 (respectively < t < 1). 
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The assumption (|2,6p is a very weak assumption on the asymptotic behavior of (j) at infinity. All 
the examples (in (i), (iii) and (v), we need to take a < 2) above Lemma [2J] satisfy this assumption. 
In fact they satisfy the following stronger assumption 

( f>(\) = \^e(X), (2.7) 

where £ is a function slowly varying at infinity. By inspecting the table in |34[ Chapter 15], one 
can come up with a lot more examples of complete Bernstein functions satisfying this stronger 
assumption. In the next example we construct a complete Bernstein function satisfying (|2.6p . but 
not the stronger (|2.7|) . 

Example 2.8 Suppose that a G (0,2). Let F be a function on [0, oo) defined by F(x) = on 
< x < 1 and 

F(x) = 2 n , 2 2 ( n ~ 1 )/ Q < x < 2 2n / a , n = 1, 2, . . . . 
Then clearly F is non-decreasing and x a / 2 < F(x) < 2x a l 2 for all x > 1. This implies that for all 

^<Hminf^M<limsu P ^M<2^. 
2 z->oo F(x) 2-i.oo F(x) 

If F were regularly varying, the above inequality would imply that the index was a/2, hence the 
limit of F{tx)/F{x) as x — > oo would be equal to ct a / 2 for some positive constant c. But this 
does not happen because of the following. Take t = 2 2// ° and a subsequence x n = 2 2n l a . Then 
tx n = 2 2 ( n+1 )/ a and therefore 

F{tx n )/F{x n ) = 2 n+2 /2 n+1 = 2 

which should be equal to ct a l 2 = c(2 2 l a ) a l 2 = 2c, implying c = 1. On the other hand, take 
any t G (1,2 2 / Q ) and the same subsequence x n = 2 2n / a . Then tx n G [2 2n / a , 2 2 (™+ 1 )A*) implying 
F{tx n ) = F(x n ). Thus the quotient F(tx n )/F(x n ) = 1 which should be equal to ct a = t a for all 
t G (1, 2 1 /"). Clearly this is impossible, so F is not regularly varying. This also shows that F(x) is 
not ~ to any cx a / 2 , as x — > oo. 

Let a be the measure corresponding to the nondecreasing function F (in the sense that a(dt) = 
F(dt)): 

oo 

o := S y^ j 2 n 5 2 2n/a ■ 
11=1 

Since Jj ^(1 + t) a(dt) < oo, a is a Stieltjes measure. Let 

2 n 



KA):=/ T T7 ^(*) = E 
■/(0,oo) A + * t=t 



X + 2 2n / a 



be the corresponding Stieltjes function. It follows from [3, Theorem 1.7.4] or [451 Lemma 6.2] that 
g is not regularly varying at infinity. Moreover, since F(x) x x a ^ 2 , x — > oo, it follows from |45[ 
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Lemma 6.3] that g(X) X A Q / 2 1 , A — > oo. Therefore, the function /(A) := l/g(X) is a complete 
Bernstein function which is not regularly varying at infinity, but satisfies /(A) x A 1_Q / 2 , A — )• oo. 

Now we are going to establish the asymptotic behaviors of u and \i under the assumption (H). 
First we claim that under the assumption (|2.6p . there exist 5 £ (0, 1) and a, so > such that 

<j)(\t) > a\ 5 (p(t), \>l,t>l/s . (2.8) 

Indeed, by Potter's theorem (cf. [31 Theorem 1.5.6]), for < e < a/2 there exists t\ such that 

«*> <2n^((±)',(^Y\=2X<, A> M > fl . 



£(At) ~ \\\tj '\t 

Hence, 

0(Ai) > c 2 {Xt) a/2 £(Xt) = c 2 t a/2 £{t)X a/2i -^- > c^{t)X a/2 -" , A > l,t > t x . 
By taking 5 := a/2 — e G (0, 1), a = C3, and sq = V*i we arrive at (|2.8p . 

Theorem 2.9 Zei S be a complete (possibly killed) subordinator with Laplace exponent <p satisfying 
(H) . TTien i/ie potential density u of S satisfies 

^xrVr 1 )- 1 ^, i^0+. (2.9) 



Proof. Put 



u;(t) : = / e- st u(s)ds ' 



Hty 

then by f)2.8|) we have 

w(Xt) < a^X^wit), A>l,i>l/s . 
Applying the second part of Proposition 12.51 we see that there is a constant c > such that 

u(t) > ci -1 ^ -1 ), 

for small t > 0. Combining this inequality with (|2.5p we arrive at (|2.9p , □ 

Theorem 2.10 Zei S be a complete subordinator with Laplace exponent (f> with zero killing coeffi- 
cient satisfying (H) . Then the Levy density [i of S satisfies 

^xrVtr^r/ 2 - 1 ^- 1 ), t^o+. (2.10) 
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Proof. Since is a complete Bernstein function, the function ip(X) := X/(p{X) is also a complete 
Bernstein function and satisfies 

A l-a/2 

where a £ (0,2) and £ are the same as in (|2.6[) . It follows from Corollary 12.41 that the potential 
measure of a killed subordinator with Laplace exponent ip has a complete monotone density v given 
by 

/oo 
[i(s)ds. 

Applying Theorem 12.91 to ip and v we get 

fj,(t,oo)=v{t)xr 1 tp(t~ 1 )~ 1 = (p(t~ 1 ), t->0. (2.11) 

By using the elementary inequality 1 — e~ cy < c(l — e~ y ) valid for all c > 1 and all y > 0, we get 
that 0(cA) < c^(A) for all c > 1 and all A > 0. Hence ^>(s _1 ) = </>(2(2s) _1 ) < 2^((2s)" 1 ) for all 
s > 0. Therefore, by (j2TTj) . for all s G (0, 1/2) 

v{s) < ci0(s _1 ) < 2ci^((2s) _1 ) < c 2 t>(2s) 

for some constants c±,C2 > 0. Since 

u(t/2) > v(t/2) - v(t) = [ n(s)ds > (t/2)fj,(t) , 

Jt/2 

we have for all t € (0, 1), 

< 2t~ 1 v(t/2) < c 2 t~ 1 v(t) < c^ 1 ^ 1 ) , 

for some constant C3 > 0. 

Using (|2.8p we get that for every A > 1 

^(s- 1 ) = 0(A(As)- 1 ) > aA^((As)" 1 ) , s < ^ . 

A 

It follows from (|2.1ip that there exists a constant C4 > such that 

C4 1 <P{s~ 1 ) < V(s) < C40(s _1 ) , S < 1 . 

Fix A := 2 1 /' 5 (( c 2 a -i) V I) 1 / 5 > 1. Then for s < (s A 1)/A, 

u(As) < C4<A((As)~ 1 ) < c 4 a- 1 A-' 5 (/)(s- 1 ) < cja^X^vis) < -v(s) 



by our choice of A. Further, 



Xs 1 , , 1 



(A - l)sfi(s) > / n{t) dt = v(s) - v(Xs) > v(s) - -v(s) = -v(s) . 

J s 2 2 

This implies that for all small t 

M(t) > 2 (\-l) t ~ 1 ' U ^ = C5t ~ lv ^ ~ c 6^V(i _1 ) 
for some constants C5, cq > 0. The proof is now complete. □ 
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3 Subordinate Brownian motion 



3.1 Definitions and technical lemma 

Let B = (Bt,F x ) be a Brownian motion in M. d with transition density p(t,x,y) = p(t,y — x) given 
by 

/ ..12 \ 

id 



p{t, x) = (ATTt)- d / 2 exp » t > 0, x,y £ 



The semigroup (P t : t > 0) of B is defined by P t f(x) = E x [f(B t )] = f Rd p(t, x, y)f(y) dy, where / is 
a nonnegative Borel function on M d . Recall that if d > 3, the Green function G^ (x, y) = G^ (x—y), 
x, y £ M. d , of B is well defined and is equal to 



G {2 \x) = / p(t,x)dt 



r(d/2-i) d+2 

Let 5 = (St : i > 0) be a complete subordinator independent of B, with Laplace exponent 
4>(X), Levy measure ^ and potential measure U. In the rest of the paper, we will always assume 
that S is a complete subordinator whose killing coefficient is zero, is dependent of B and satisfies 
(H). Hence lim^oo 0(A) = oo, and thus S has a completely monotone potential density u. We 
define a new process X = (Xt : t > 0) by Xt := -Bs*- Then X is a Levy process with characteristic 
exponent <&(x) = 0(|x| 2 ) (see e.g.[33j pp. 197-198]) called a subordinate Brownian motion. The 
semigroup (Qt : i > 0) of the process X is given by 

POD 

Q t f{x) = E x [f{X t )] = E x [f(B St )\ = / P s f(x) P(S t e ds) . 

Jo 

The semigroup Qt has a density q(t, x, y) = q(t, x — y) given by q(t, x) = f Q °° p(s, x) F(St € ds). 

Recall that, according to the criterion of Chung-Fuchs type (see [30] or [33], pp. 252-253]), X 
is transient if and only if for some small r > 0, fi x i <r dx < oo. Since = </>(|2;| 2 ), it follows 
that X is transient if and only if 

dA<oo. (3.1) 



<P(X) 

This is always true if d > 3, and, depending on the subordinator, may be true for d = 1 or d = 2. 
In the case d < 2, if there exists 7 E [0, d/2) such that 

liminf^>0, (3.2) 

then ([3T]) holds. 

For x G M d and a Borel subset j4 of R , the potential measure of X is given by 

roo roo /*oo poo 

G(x,A) = E x l {XteA} dt= Q t l A (x)dt= / P s lA(x)F(S t £ds)dt 
Jo Jo Jo Jo 

f'OG p roo 

= / P s lAu(s)ds= / p(s,x,y)u(s)dsdy , 

io JA Jo 
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where the second line follows from (|2.ip . If X is transient and A is bounded, then G(x,A) < oo 
for every x S R . In this case we denote by G{x,y) the density of the potential measure G(x,-). 
Clearly, G(x, y) = G{y — x) where 

/'OO poo 

G{x)= I p(t,x)U(dt) = / p{t,x)u(t)dt. (3.3) 
Jo Jo 

The Levy measure II of X is given by (see e.g. [331 PP- 197-198]) 

p poo p 

IL(A) = pit, x) n(dt) dx = J(x) dx, 4cR d , 

J A Jo J A 

where 

f oo poo 



J{x):= I p{t,x) n{dt) = I p(t,x)fi(t)dt (3.4) 
Jo Jo 

is the Levy density of X. Define the function j : (0, oo) — > (0, oo) by 

i (r) : =^°°(4vr)- d / 2 t- d / 2 e X p(-^ »(dt) , r>0, (3.5) 

and note that by ([33D, J(x) = j(\x\), x £ R d \ {0}. 

Since x i-> p(t, x) is continuous and radially decreasing, we conclude that both G and J are 
continuous on R d \ {0} and radially decreasing. 

The following technical lemma will play a key role in establishing the asymptotic behaviors of 
the Green function G and the Levy density J of the subordinate Brownian motion X in the next 
subsection. 

Lemma 3.1 Suppose that w : (0, oo) — > (0, oo) is a decreasing function, £ : (0, oo) — > (0,oo) a 
measurable function which is locally bounded above and below by positive constants and is slowly 
varying at oo, and (3 & [0, 2], /3 > 1 — d/2. If d = 1 or d = 2, we additionally assume that there 
exist constants c > and 7 < d/2 such that 

w{t)<ctf- 1 , Vt>l. (3.6) 

Let 

roo 1 ,2 

I(x):= I {4irt)- d / 2 e-^Tw(t)dt. 

N If 







then 

1 w(\x\ 2 ) 
I t) x x — r — > f) 

I I \\x\ 2 / 1 1 
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(b) If 

<Q~m/r)' (3 - 8) 

then 

, . T(d/2 + B-l) 1 

J W 4 l-/V*/2 1x1^-2^^)' R-^ U ' 

Proof, (a) Let us first note that the assumptions of the lemma guarantee that I(x) < 00 for every 
x 7^ 0. Now, let £ > 1/4 to be chosen later. By a change of variable we get 



(W*'f,W + N-' + %w) ■ (3.9) 



We first consider Ii(x) for the case d = 1 or d = 2. It follows from the assumptions that there 
exists a positive constant c\ such that w(s) < cis 7_1 for all s > l/(4£). Thus 

/,(,•)£ /^'V^e-'dfSy di<c 2 |x| 2 ^ 2 / ?N i^-T- 1 eft = 031x1^- 



It follows that 



Mir^M^rf-'^))^. (3,0) 



In the case <i > 3, we proceed similarly, using the bound w;(s) < u;(l/(4£)) for s > l/(4£). 
Now we consider hi^x): 

|2- 



ixr"+ 2 



Using the assumption (|3.7p . we can see that there is a constant c\ > 1 such that 

for all t and x satisfying |x| 2 /(4t) < l/(4£). 

Now choose a 5 G (0, d/2 — 1 + /3) (note that by assumption, d/2 — 1 + /3 > 0). By Potter's 
theorem (cf. [31 Theorem 1.5.6 (i)]), there exists p = p{5) > 1 such that 

^ <2 ((V\*\*\\.(V\*\*Y* 



■ At 



(iff) v (WRF) j = 2 (< 4 " 4 v W) < ^ v r°) (3,1) 
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whenever ^jp > p and j^p > p. By reversing the roles of l/|x| 2 and 4i/|x| 2 we also get that 



for T^p > p and A > p. Now we define £ := | so that for all x 7^ with \x\ 2 < ^ and t > 
we have that 



|x| 

< Cl c 2 t d / 2 - 2+ Pe-\t 5 vr 5 ). (3.13) 



Let 

/>oo 

c 3 := q 1 ^ 1 / t d / 2 - 2 ^e-\t 8 A r 5 )di < 00 , 
Jo 

/•oo 

c 4 := cic 2 / i d / 2 ~ 2+/ V(^ V^ 5 )(ft < 00. 
Jo 

The integrals are finite because of assumption d/2 — 2 + f3 — 5 > — 1. It follows from (|3.13|) that 



|z|->0 JO 

This means that 



- d + 2 h(x) (\x\ d - 2r >+ 2 t(jL)\ 



* I - ( !£) (!£)' t («, ) g> * x !. ,3,4, 



Combining (|3.1U|) and (|3.14|) we have proved the first part of the lemma. 

(b) The proof is almost the same with a small difference at the very end. Since t is slowly varying 
at 00, we have that 

lim^ = l. 

This implies that 

|2\ /l™!2\/3 ' • - 1 
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By the right-hand side inequality in (|3.13p . we can apply the dominated convergence theorem to 
conclude that 

lim \x\- d+2 I 2 {x) ( \ x \ d -^+H{-^)\ 
W->o V \ x \ J 

= 4 /3 r(d/2 - 1 + p) . 

Together with (13, and A3. 10f> this proves the second part of the lemma. □ 



3.2 Asymptotic behavior of the Green function and Levy density 

The goal of this subsection is to establish the asymptotic behaviors of the Green function G(x) 
and Levy density J(x) of the subordinate process X under certain assumptions on the Laplace 
exponent <p of the subordinator S. We start with the Green function. 

Theorem 3.2 Suppose that the Laplace exponent <p is a complete Bernstein function satisfying the 
assumption (H) and that a £ (0, 2 A d). In the case d < 2, we further assume (|3.2p . Then 



^ \x\ d (j>(\x\- 2 ) ^ \x\ d - a l{\x\-^ |j | ' 



Proof. It follows from Theorem 12.91 that the potential density u of S satisfies 

£(t- 

Using (|2.5|) and (|3.2|) we conclude that in case d < 2 there exists c > such that 



U (t) x t- i <p{t- i r i x — — , t^o+. 



u(i)<cf 7-1 , t>l. 

We can now apply Lemma 13.11 with w(t) = u(t), (3 = 1 — a/2 to obtain the required asymptotic 
behavior. □ 



Remark 3.3 (i) Since a is always assumed to be in (0,2), the assumption a £ (0,2 Ad) in the 
theorem above makes a difference only in the case d = 1. 

(ii) In case d > 3, the conclusion of the theorem above is proved in [46\ Theorem 1 (ii)-(iii)] under 
weaker assumptions. The statement of |46l Theorem 1 (ii)] in case d < 2 is incorrect and the proof 
has an error. 

The asymptotic behavior near the origin of J(x) is contained in the following result. 
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Theorem 3.4 Suppose that the Laplace exponent cj) is a complete Bernstein function satisfying the 
assumption (H). Then 



Proof. It follows from Theorem 12. 101 that the Levy density /i of S satisfies 

n(t) x x r / 2 - 1 ^- 1 ) , f -> o + . 

Since fi(t) is decreasing and integrable at infinity, one can easily show that there exists c > such 
that 

fi(t) < cr 1 , t > l. 

In fact, if the claim above were not valid, we could find an increasing sequence {t n } such that 
t\ > 1, t n t oo, t n — t n _i > t n /2 and that n(t n ) > nt' 1 . Then we would have 

roo rti 00 rt n fl — 1 °° jj 

/ M (t)«ft = / n(t)dt + v / iM(t)dt > ^— +y,^=oo, 

h h „=2- /t «-i *! n=2 2 

contradicting the integrability of /i at infinity. Therefore the claim above is valid. We can now 
apply Lemma 13.11 with w(t) = fi(t), f3 = 1 + a/2 and 7 = to obtain the required asymptotic 
behavior. □ 



Proposition 3.5 Suppose that the Laplace exponent (f> is a complete Bernstein function satisfying 
the assumption (H) . Then the following assertions hold. 

(a) For any K > 0, there exists C4 = C^{K) > 1 such that 

j(r) <C 4 j(2r), Vre(0,tf). (3.15) 

(b) There exists C5 > 1 such that 

j(r)<C 5 j(r + l), Vr>l. (3.16) 

Proof. (|3.15p follows immediately from Theorem 13.41 However, we give below a proof of both 
(|3~T5jl and (l3TT6|) using only (l3TT7D - (^T8]l . 

For simplicity we redefine in this proof the function j by dropping the factor (47r) — d / 2 from its 
definition. This does not effect (I3.15P and (I3.16p . It follows from Lemma 1 2. II and Theorem 12.101 
that 

(a) For any K > 0, there exists c\ = c\ (K) > 1 such that 

fj,(r) < cxfj,(2r), Vr e (Q,K). (3.17) 
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(b) There exists C2 > 1 such that 

H{r) < C2/i(r + l), Vr>l. (3.18) 

Let < r < K. We have 

/»oo 

j(2r)= t- d / 2 exp(-r 2 /t)n(t)dt 



> \(^j Q ° t~ d/2 exp(-r 2 /t)//(t) dt + jf** r d / 2 exp(-r 2 /t)fi(t) (ttj 



\{h+h)- 



Now, 

poo 2 /-co 2 o 2 

h = / t- d / 2 eM--)Kt)dt= i - d / 2 exp(--)exp(-— )//(*) eft 

/•oo 2 o r 2 /-oo 2 

> / t^/ 2 exp(- )exp(-i-V(0^> e - 3 ^ 2 / t - d / 2 eM- T Mt)dt. 

JK/2 4t 7iC/2 4 * 

h = t- d / 2 exp(- — )^(t)dt = 4~ d / 2+1 / s- d/2 exp(-^-M4s)ds 

Jo * Jo 4s 



-2 4 -d/2+l / -d/2, 



> c{ 2 A- a l^ i s - d /^ exp( __ )/i(s)ds . 



Combining the three displays above we get that j(2r) > csj(r) for all r 6 (0, if). 
To prove (|3.16p we first note that for all t > 2 and all r > 1 it holds that 

(r + 1) 2 r 2 

- — ! — < 1 . 

t t-1 ~ 

This implies that 

( r + l)2\ / r 2 



exp 
Now we have 



4/ 



} > e~ 1/4 exp (- 4 ^ fj ) , for all r > l,t > 2. (3.19) 



j(r + 1) = t^/ 2 exp(-^±^)/i(t) (it 



> 
~ 2 



I t- d ' 2 «p(-£±^-)/x(t) + t^/ 2 exp(-^-) M (t) (it 



i(^3 + / 4 ). 
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For I 3 note that (r + l) 2 < 4r 2 for all r > 1. Thus 

'8 ^in2 



= / t~ d/2 exp(- ^ + ^ )^{t)dt > [ t- d / 2 eM-r 2 /t)fi(t)dt 
Jo 4t Jo 

f 2 r 2 
= 4- d / 2+1 y s - d / 2 exp(-— )/i(4s)ds 

> cr 2 4- d / 2+1 jT s - rf / 2 exp(-^)/x( S ) , 



> / t- d/2 exp{-l/4)exp(-— — ) u(t) dt 



3 



/ r^exiK-l/^expC-^-^)^). 

/•oo 2 

= e" 1/4 y (s-l)- d/2 exp(-— )ii(s + l)ds 

r°° r 2 
> c^e- l/4 J s- d/2 exp(-—)fi(s)ds. 

Combining the three displays above we get that j(r + 1) > C4 j{r) for all r > 1. □ 
3.3 Some results on subordinate Brownian motion in R 

In this subsection we assume d = 1. We will consider subordinate Brownian motions in R. Let 
B = (Bt : t > 0) be a Brownian motion in R, independent of S", with 



E 



e iO(B t -B ) 



The subordinate Brownian motion X = (X : i > 0) in R defined by X^ = Bg t is a symmetric 
Levy process with the characteristic exponent = (ft(9 2 ) for all 6 G R. In the first part of this 
subsection, up to Corollary 13.81 we do not need to assume that <\> satisfies the assumption (H). 

Let X t := sup{0 V X s : < s < t} and let L = (L t : t> 0) be a local time of X - X at 0. L 
is also called a local time of the process X reflected at the supremum. Then the right continuous 
inverse L^ 1 of L is a subordinator and is called the ladder time process of X. The process X r -i 
is also a subordinator and is called the ladder height process of X. (For the basic properties of 
the ladder time and ladder height processes, we refer our readers to [2, Chapter 6].) Let x De t ne 
Laplace exponent of the ladder height process of X. It follows from [201 Corollary 9.7] that 

The next result, first proved independently in [27\ and [28J, tells us that x 1S a ^ so complete Bernstein 
function. The proof presented below is taken from [27] . 
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Proposition 3.6 Suppose cf>, the Laplace exponent of the subordinator S, is a complete Bernstein 
function. Then the Laplace exponent x of the ladder height process of the subordinate Brownian 
motion X t = B$ t is also a complete Bernstein function. 

Proof. It follows from Theorem |34} Theorem 6.10] that <f> has the following representation: 

1 „ g «A) = 7 + ^( T ^- x ^),( t ) <i( , (3.21) 

where r\ in a function such that < n(t) < 1 for all t > 0. By (|3.2ip and (|3.20p . we have 

7 1 f°° f°° ( t 1 \ . w d9 

log X (A) = I + -/ / l-—-——)r,(t)dt. 



2 nJo J \i + t 2 \ 2 e 2 + t) IK> i + e 2 

By using < r/(t) < 1, we have 
t 1 



r/(i) 



i + 1 2 \ 2 e 2 + 1 



i 1/1 . \ 2 e 2 t 



i + e 2 ~ i + t 2 i + o 2 \\ 2 e 2 + t \ 2 o 2 + t 
1/1 \ 2 t 

< z : To . . + 



i + t 2 \ x 2 e 2 + i A 2 <9 2 + 1 



Since 



1 ,„ 1 f°° 1 ,„ 1 y/t 



d0=~ T2a2 7 d6 = — — — / d 7 



A 2 2 + f ^ + 1 t A Jo 7 2 + 1 2A^t 



7T 



we can use Fubini's theorem to get 



logx(A) - l + f(^L_-^-^), m (3.22) 

s 1 
1 + s 2 ~ A + s 



71 + / ( ^2 - — ) V(s 2 )ds. 



Applying [341 Theorem 6.10] we get that x 1S a complete Bernstein function. □ 

The potential measure of the ladder height process of X is denoted by V and its density by v. 
We will also use V to denote the renewal function of X: V(t) := V((0, t)) = Jq v(s) ds. 
The following result is first proved in |27j . 



Proposition 3.7 x ^ s related to <fi by the following relation 

e~ V V<K^ 2 ) < X(A) < e" /2 y/</>(\ 2 ) , for all A > . 
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Proof. According to (|3.22p . we have 

Together with representation (|3.2ip we get that for all A > 



rj(t)dt. 



lo gx (A)--log0(A 2 ) 



f((r 



< 



1 

2 ,„ 



+ t 2 ^i(\ + ^t)J \l + t 2 \ 2 + t 

A( v / t + A) 1 f°° A 

■dt = - — = dt 



rj(t) dt 



2 J (A 2 + t)y/i(X + 2j (\ 2 + t)yft 



IT 



This implies that 



i.e., 



-vr/2 < logx(A) - -log0(A 2 ) < vr/2, for all A > , 



-* /2 < x w<k* 



2^-1/2 < e n/2 



for all A > . 



□ 



Combining the above two propositions with Corollary 12.31 we obtain 

Corollary 3.8 Suppose <f>, the Laplace exponent of the subordinator S, is a complete Bernstein 
function satisfying lim^oo 0(A) = oo. Then the potential measure of the ladder height process of 
the subordinate Brownian motion Xt = B$ t has a completely monotone density v. In particular, v 
and the renewal function V are C°° functions. 

In the remainder of this paper we will always assume that <p satisfies the assumption (H). We 
will not explicitly mention this assumption anymore. 

Since <f>(\) X A Q / 2 ^(A) as A — > oo, Lemma [3771 implies that 

X (\)~\ a/2 (l(\ 2 )) 1/2 , t^oo. (3.23) 

It follows from (13.23j) that lim^oo x(A)/A = 0, hence the ladder height process does not have a 
drift. Recall that V(t) = V((0,t)) = J^v(s)ds is the renewal function of the ladder height process 



of X. In light of ()3.23p . we have, as a consequence of Theorem 
Proposition 3.9 As t — > 0, we have 

v(t) x ^(t- 2 )- 1 /' 2 

and 



the following result. 



t 



a/2 



v(t) x rVr 2 )- 1/2 



(e(t' 2 )y/ 2 

t a ' 2 - 1 
{i{t- 2 )f/ 2 
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Remark 3.10 It follows immediately from the proposition above that there exists a positive con- 
stant c> such that V(2t) < cV(t) for all t G (0, 2). 

It follows from (|3.23[) above and [29\ Lemma 7.10] that the process X does not creep upwards. 
Since X is symmetric, we know that X also does not creep downwards. Thus if, for any o 6 R, we 
define 

T a = inf{t > : X t < a}, a a = inf{t > : X t < a}, 

then we have 

Pz(r a = oa) = 1, x > a. (3.24) 
Let Gr 0iOO \(x,y) be the Green function of X in (0,oo). Then we have the following result. 

Proposition 3.11 For any x,y > we have 



G 



(0,oo) 



Jq v(z)v(y + z - x)dz, x < y, 
fx-y v ( z ) v (y + z- x )dz, x>y. 



Proof. Let _XT<o,oo) be the 

process obtained by killing X upon exiting from (0, oo). By using (|3.24[) 
above and [21 Theorem 20, p. 176] we get that for any nonnegative function on / on (0, oo), 



E, 



roc 

/ f(xi°^)dt 

Jo 



^0 



v(z)f(x + z - y)v(y)dzdy , 



(3.25) 



where k is the constant depending on the normalization of the local time of the process X reflected 
at its supremum. We choose k = 1. Then 



E, 



Jo 



v (y) v(z)f(x + y- z)dzdy 



v{z) / v(w + z — x)f(w)dwdz 



v ( z ) / v(y)f(x + y - z)dydz 



f{w) I v(z)v(w + z — x)dzdw + I f(w) / v{z)v{w + z — x)dzdw 



On the other hand, 



E, 



roo 

/ f{X?>^)dt 
Jo 



G^oo)(x,w)f(w)dw. 



(3.26) 



(3.27) 



□ 



By comparing (|3.26p and (|3.27p we arrive at our desired conclusion. 

For any r > 0, let Gr(o,r) be the Green function of X in (0, r). Then we have the following result 
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Proposition 3.12 For all r > and all x G (0, r) 

G {Q>r) {x,y)dy<2V{x)V{r). 



In particular, for any R > 0, there exists Cq = Cq(R) > such that for all r G (0,i?) and all 
x G (0,r), 



r «/2 x a//2 



G(o, r )(x,y)dy < C 6 (0(r 2 )#E 2 )) V2 X (£(r - 2)) l/2 (l (a .-2))-l/2 

Proof. For any a; G (0, r), we have 



G(o,r){x,y)dy < / G^ oo) (x,y)dy 
o ■/ o 



v(z)v(y + z — x)dzdy + / / v(z)v(y + z — x)dzdy 
o ix Jo 

u(z) / v(y + z — x)dydz + / w(z) / v(y + z — x)dydz < 2V(r)V(x). 
o A-z Jo A 

Now the desired conclusion follows easily from Proposition 13.91 □ 

As a consequence of the result above, we immediately get the following. 
Corollary 3.13 For all r > and all x G (0, r) 

G (0 , r) (x, y) rig < 2F(r) (V(x) A V(r - x)) . 
In particular, for any R > 0, i/iere exists Cj = Cj(R) > suc/i i/iai /or a// x G (0, r), and r G (0, R), 







G { v tT) {x,v)dy < C 7 (0(r- 2 ))- 1 /2^ (0(:r ^ )r i/2 A((/)((r ._ 3;r 2 )r i/2^ 



r a/2 / j.a/2 /_ _ x )a/2 

I A 



(£(r- 2 ))V2 ^(^( x -2))l/2 (^((r-aT)-2))i/2y ' 

Proof. The first inequality is a consequence of the identity J" Q r Gr 0r \(x, y)dy = L G( 0jr )( r — x -, V)dy 
which is true by symmetry of the process X. The second one now follows exactly as in the proof 
of Proposition 13.121 □ 



Remark 3.14 With self-explanatory notation, an immediate consequence of the above corollary 
is the following estimate 

•r 

G(-r,r) (»> V) d y < 2V(2r) [V(r + x) A V{r - x)) . (3.28) 
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4 Harnack inequality and Boundary Harnack principle 



From now on we will always assume that X is a subordinate Brownian motion in M. d . Recall that 
(H) is the standing assumptions on the Laplace exponent <p. The goal of this section is to show 
that the Harnack inequality and the boundary Harnack principle hold for X. The infinitesimal 
generator L of the corresponding semigroup is given by 



for / G C 2 



L/(x) = / (f(x + y) - f{x) - y ■ V/(x)l { |„|< 1} ) Jiy)dy 
Moreover, for every / G C 2 (M d ) 

ft 

\ ds 



(4.1) 



f(X t ) - f(X ) - f\f(X s 
Jo 



is a P x -martingale for every x G M. d . We recall the Levy system formula for X which describes 
the jumps of the process X: for any non- negative measurable function / on R + x R rf x W 1 with 
f{s,y,y) = for all y G M. d , any stopping time T (with respect to the filtration of X) and any 

x G M. d , 

(4.2) 



E, 



Y,fi*,Xs-,X s 



s<T 



E, 



J Uf(s,X s ,y)J(X s ,y)dyjds 



(See, for example, \15\ Proof of Lemma 4.7] and \16\ Appendix A].) 
4.1 Harnack inequality 

It follows from Theorem 13.41 and the 0- version of [3J Propositions 1.5.8 and 1.5.10] that 



and 



r- 2 / s d+1 j(s)ds x ^ 
Jo r 



°° s d -'j{s)ds^ 1 -^ 



r -»■ 0. 



(4.3) 



(4.4) 



For any open set D, we use td to denote the first exit time from D, i.e., td = inf{t > : X< ^ 

D}. 

Lemma 4.1 There exists a constant C$ > swc/i t/iat /or every r G (0, 1) and every t > 0, 

x ^sup |J*G - X | > rj < Cs<Pir~ 2 )t . 

Proof. It suffices to prove the lemma for x = 0. Let / G C fe 2 (M d ), < / < 1, /(0) = 0, and /(y) = 1 
for all \y\ > 1. Let Cl = sup y E j)fe l(^ 2 /%^)/(2/)l- Then |/(z + y) - f(z) - y ■ Vf(z)\ < f \y\ 2 . 
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For r € (0, 1), let f r (y) = f{y/f). Then the following estimate is valid: 

ci \y? 

\f r (z + y)- f r (z) - y ■ V f r (z)l { \y\< r} \ < — — l{|„|< r } + 1 {\y\>r} 

< c 2(l{| J/ |<r}^i- + l{|y|>r})- 

By using (|4.3p and (|4.4p . we get the following estimate: 

\lfr(z)\ < [ \fr(z + y)-fr(z)-yVfr(z)l M <r)\J(y)dy 

< c 2 f \Hy\<r}-Za +1 {\y\>r}] J (y)dy 

< C 8 0(r- 2 ), (4.5) 
where the constant Cg is independent of r. Further, by the martingale property, 

Eo/r(X Tfl(0ir)At ) - / r (0) = E / Lf r (X s ) ds (4.6) 

J 

implying the estimate 

E f r (X TB(O r)M ) < CsHr' 2 )t . 

If X exits 5(0, r) before time t, then f r (X TB(0 r)At ) = 1, so the left hand side is larger than 
n(r B (o,r) <t). □ 

Lemma 4.2 For every r £ (0, 1), anc? every x 6 M d , 

inf E^ r B(a . )] > - , 
where C 8 is the constant from Lemma \4-l\ 

Proof. Using f|4.5j) and (|4.6p we get that for any t > and z 6 B(x,r/2), 

Mr B (o,r/2) <t) < C 8 0((r/2)^ 2 )E O [r B(0)r/2) A t] 
= C 8 0((r/2)- 2 )E 2 [re^/z) A t] 
< C 8 0((r/2)- 2 )E, [r B(:Cjr) A t] . 

Letting i —> oo, we immediately get the desired conclusion. □ 

Lemma 4.3 There exists a constant Cg > such that for every r E (0, 1) and every x £ M d , 

sup E, [r B( , ir) ] < -^K-y 
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Proof. Let r S (0, 1), and let x S M. d . Using the Levy system formula (|4.2p . we get 

1 > F z (\X TB{x r) - x\ > r) 



G B (x,r)(z,y) / j(\u-y\)dudy , 

B{x,r) J B(x,r) 

where G B ( x ,r) denotes the Green function of the process X in B(x,r). Now we estimate the inner 

C 

integral. Let y E B(x, r), u G B(x, r) . If u € B(x, 2), then |« — y\ < 2\u — x\, while for u ^ B(x, 2) 
we use | it — y\ < \u — x\ + 1. Then 



/ J(W-y\)du 

JB(x,r) 



B(x,r) nB(a;,2) JB(x,r) C\B{x,2) 



j(\u-y\)du+ / j(\u-y\)du 



> / j(2|it-a;|)du+ / j(\u - x\ + 1) du 

J B(x,rfr\B(x,2) J B(x,r) C nB(x,2) c 

> / c~~ j(\u — x\) du + / c~ j(\u — x\) du 

J B(x,rfnB(x,2) J B(x,r) C nB(x,2) c 

c j(\u — x\) du , 



<B{x,r) 

where in the next to last line we used ()3. 15[) and (|3.16p . Now, It follows from (|4.4p that 

1 > / G B{x ^ r) (z,y)dy / c~ l j{\u - x\) du 

JB{x,r) JB{x,r) 
roo 

= E z [T B (x,r)] C _1 Ci / (fo 
ir 

= c 2 4>(r' 2 )E z [r B ( x ,r)] 

which implies the lemma. □ 

An improved version of the above lemma will be given in Proposition 14.91 later on. 

Lemma 4.4 There exists a constant C\$ > such that for every r £ (0, 1), every x G M. d , and any 
A C B(x,r) 

\A\ 

F y (T A < T B(X:3r) ) > C 10 — L-!— , for allyeB(x, 2r) . 

Proof. Without loss of generality assume that ¥ v {Ta < T B ( x ^ r \) < 1/4. Set r = T B ( x ^ r y By 
LemmaSH P y (r < t) < P y (T B(y ^ < t) < ci<j){r- 2 )t. Choose t = 1 / ' (Ac^r- 2 )) , so that Pj,(t < 
*o) < 1/4- Further, if z G -B(x,3r) and m £ A C B(x,r), then |u — z| < 4r. Since j is decreasing, 
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j(\u — z\) > j(4r). Thus, 



,(Ta<t) > E y 1 {x s ^x s ,x s eA} 

s<T a AtM 
rT A At At f 

= E y / / j{\u- X s \)duds 



J A 

T A f\T?\t Q 



> Ey / j{Ar)duds 

JO J A 

= j(4r)\A\E y [T A ATAto}, 



where in the second line we used properties of the Levy system. Next, 

E v [Ta A t A to] > E y [t ;T A >T>t ] 

= t Fy(T A >T>t ) 

> t [l-F y (T A <T)-Fy(T<t )] 

> tl 1 



2 8c 1( j)(r- 2 ) 



■i 

The last two displays give that 



1 1 , „ 3(4r) 



\(TA<r)> j (Ar)\A\^^ ) =-\A\^ ) 



The claim now follows immediately from (|2.6p and Theorem 13.41 □ 

Lemma 4.5 There exist positive constant C\\ and C\2, such that if r G (0, 1), x £ M. d , z £ B(x,r), 
and H is a bounded nonnegative function with support in B(x,2r) c , then 

E z H(X TB ^ r) )<C n E z [T B{Xtr) ] J H(y)j(\y-x\)dy, 

and 

E z H(X TB{x r) )>C 12 E z [T BM \ j H(y)j(\y-x\)dy. 

Proof. Let y € B(x, r) and u G B(x, 2r) c . If u £ B(x, 2) we use the estimates 

2~ 1 |'ii - x\ < \u - y\ < 2\u - x\, (4.7) 

while if u ^ B(x, 2) we use 

\u — x\ — 1 < \u — y\ < \u — x\ + 1. (4-8) 
Let .B C 2r) c . Then using the Levy system we get 



E, 



1b(X t ) =E Z / j(\u-X s \)duds. 
' J Jo Jb 
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By use of (|3.15[) . (|3.16p . (|4.7p . and (|4.8]) . the inner integral is estimated as follows: 



j(\u-X s \)du = / j(\ u - Xs \)du+ / j{\u-X s \)du 

B JBnB(x,2) JBnB(x,2) c 



< / j{2' 1 \u- x\)du + I 

JBnB(x,2) JE 



j(\u — x\ — 1) dii 



BnB(x,2) c 

< / cj{\u — x\) du + / cj(|u — x|)cfai 

iflnB(i,2) iflnB(i,2) c 

= c / — x|) (in. 



Therefore 



E, 



1b (X 



T B(x,r) > 



< E, 



T B(x,r) 



c I j(\u — x\) du 

' B 



C^z(T B (x,r)) J lB(u)j(\u - x\)du. 



Using linearity we get the above inequality when Is is replaced by a simple function. Approximating 
if by simple functions and taking limits we have the first inequality in the statement of the lemma. 
The second inequality is proved in the same way. □ 

Definition 4.6 Let D be an open subset o/R rf . A function u defined on M. d is said to be 

(1) harmonic in D with respect to X if 

K x [\u(X Tg )\] < oo and u(x) = K x [u(X TB )\ , x&B, 
for every open set B whose closure is a compact subset of D; 

(2) regular harmonic in D with respect to X if it is harmonic in D with respect to X and for each 

x e D, u(x) = E x [u(X TD )\ . 

Now we give the proof of Harnack inequality. The proof below is basically the proof given in 
|39j which is an adaptation of the proof given in [I]. However, the proof below corrects some typos 
in the proof given in |39j. 

Theorem 4.7 There exists C\% > such that, for any r € (0, 1/4), xo G and any function u 
which is nonnegative on M. d and harmonic with respect to X in B(xq, 17r), we have 



u(x) < Ci 3 u(y), for all x,y £ B(x ,r). 
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Proof. Without loss of generality we may assume that u is strictly positive in B(xq, 16r). Indeed, 
if u(x) = for some x G B(xq, 16r), then by harmonicity = u(x) = K x [u(X TB )] for x £ B = 
B(x,e) C B(xo,16r). This and the fact that the Levy measure of X is supported on all of M. d 
and has a density imply that u = a.e. with respect to Lebesgue measure. Moreover, by the 
harmonicity, for every y G B(xo,16r), u(y) = K y [u(X TB )] = where B = B(y,5) C B(x ,l6r). 
Therefore, if u(x) = for some x, then u is identically zero in B(xq, 16r) and there is nothing to 
prove. 

We first assume u is bounded on M. d . Using the harmonicity of u and Lemma 14.41 one can show 
that u is bounded from below on B(xq, r) by a positive number. To see this, let e > be such that 
F = {x G B(xo,3r) \ B(xo,2r) : u(x) > e} has positive Lebesgue measure. Take a compact subset 
K of F so that it has positive Lebesgue measure. Then by Lemma 14.41 f° r x £ B(xo,r), we have 



u(x) = K x u(X TkAtb 



3(x ,9r)' 



> ce 



\K\ 



B(x ,3r)\ 



(T A < r B(a . 0il6r) ) > c 2 Vy G B(x ,8r). (4.10) 



for some c > 0. By taking a constant multiple of u we may assume that xa£s(x ,r) u = 1/2- Choose 
zo G B(xo,r) such that it(^o) < 1- We want to show that u is bounded above in B(xo,r) by a 
positive constant independent of u and r G (0,1/4). We will establish this by contradiction: If 
there exists a point x G B(xq, r) with u(x) = K where K is too large, we can obtain a sequence of 
points in B(xq, 2r) along which u is unbounded. 

Using Lemmas I4.2ll4.3l and l4.51 one can see that there exists c\ > such that if x G s G (0, 1) 
and H is nonnegative bounded function with support in B(x,2s) c , then for any y,z G B(x,s/2), 

^H(X TB(x J < Cl E y H(X TB(x J. (4.9) 

By Lemma 14.41 there exists c 2 > such that if A C £?(xo,4r) then 

y[lA<r B{xo , 16r) ) > C2^- 4r)| 

Again by Lemma 14.41 there exists C3 > such that if x G s G (0, 1) and F C B(x,s/3) with 
|F|/|fi(a;,s/3)| > 1/3, then 

P. (Tp < r B{X;S) ) > ca. (4.11) 

Let 

r/=|, C = dA-)7?. (4.12) 
3 3 ci 

Now suppose there exists x G -B(xo,?") with u(x) = K for K > Kq := 2 \ B ^^ 1 ^ \j 2 ^ . Let s be 
chosen so that 

|B(x ,£)| - 3^21 < 1. (4.13) 

Note that this implies 

2 \ 1/d 



s = 12 — riT-Vd < r ^ 414 n 
c 2 C/ 
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Let us write B s for B(x, s), t s for tb( x ,s)i an d similarly for B2 S and t 2s . Let A be a compact subset 
of 

A> = {yeB(x,^):u(y)>(K}. 

It is well known that u(X t ) is right continuous in [0, T B r XOj i Sr \). Since zq G B(xq,t) and A' C 
-B(z, §) C i?(xo, 2r), we can apply (|4.1Up to get 

1 > u(zq) >^zM X T A AT B(xo , 1(ir) )l{T A <T B(xo , 16r) }} 

> (KF Z0 (T A <T B(x0!l6r) ) 
~ C2CK \B(x ,Ar)\- 

Hence 

\A\ < |B(x ,4r)| _ 1 



|fl(s, f)| - oaCK-IBCs.f)! 2 
This implies that \A'\/\B(x,s/3)\ < 1/2. Let F be a compact subset of B(x, s/3) \ A' such that 

|F| >i (4.15) 



|B(x,f)| - 3 
Let H = u ■ 1 rc . We claim that 

E x [u(X Ts );X Ts $ B 2s ]<r ] K. 

If not, E x H(X Ts ) > rjK, and by (H~9j) . for all y G s/3), we have 

u(y) = E y u(X Ta ) > E y [u(X Ts );X T3 £ B 2s ] 
> c^HiX^^c^rjK^CK, 

contradicting (I4,15h and the definition of A'. 
Let M = sup£ 2s u. We then have 

K = u(x)=E x [u(X TsATf )] 

= E x [u(X Tf );T f <t s ]+ E x [u(X Ts );t s <T f ,X Ts G B 2s ] 

+E x [u{X Ts );t s <T F ,X Ts £B 2s ] 
< (KF X (T F < t.) + M¥ x (t s < T F ) + V K 
= (KF X (T F < t 8 ) + M(l - F X (T F < t s )) + rjK, 

or equivalently 

M 1 — 7y — C > 

K " 1 - P x .(7> < r s ) + ^ 

Using (|4.1ip and (I4.12|) we see that there exists f3 > such that M > if (1 + 2/3). Therefore there 
exists x' G B(x, 2s) with u(x') > K(l + /3). 
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Now suppose there exists x\ G B(xo,r) with u(x\) = K % > Kq. Define s\ in terms of K\ 
analogously to (|4.13p . Using the above argument (with x\ replacing x and X2 replacing x'), there 
exists X2 G B(x\, 2s\) with u(x2) = K^ > (1 + 0)K\. We continue and obtain S2 and then X3, K$, 
s 3 , etc. Note that x i+1 G B^^Sj) and K { > (1 + In view of (jjHD , 

E 1^1-^1 ^ r+24 r ^ K * 

i=0 v 2W i=l 



-1/d 



f 2 \ ' 

< , + 24^-j ^ 1/rf r^(l+/?)- 

= r + 24r — K~ 1/d r ^(1 + /3)-^ 

^ C2< '^ i=0 



r + C4rA 1 



where c 4 := 24(^) 1 ^£^ (l + p)-*/*. So if i^i > V if then we have a sequence xi, X2, • • • 
contained in B(xQ,2r) with u(xi) > (1 + /3)* _1 ivTi — )■ 00, a contradiction to u being bounded. 
Therefore we can not take K\ larger than cf V Kq, and thus sup y€B ( Xo ^ u(y) < cf V Kq, which is 
what we set out to prove. 

In the case that u is unbounded, one can follow the simple limit argument in the proof of |39} 
Theorem 2.4] to finish the proof. □ 

By using the standard chain argument one can derive the following form of Harnack inequality. 

Corollary 4.8 For every a G (0, 1), there exists C14 = C\±{a) > such that for every r G (0, 1/4), 
xq G M. d , and any function u which is nonnegative on M. d and harmonic with respect to X in 
B(xq, r), we have 

u(x) < Cuu(y), for all x, y G B(xQ,ar) . 
4.2 Some estimates for the Poisson kernel 

Recall that for any open set D in M. d , td is the first exit time of X from D. 

We recall from Subsection 13.11 that X has a transition density q(t,x,y), which is jointly contin- 
uous. Using this and the strong Markov property, one can easily check that 

q D (t,x,y) := q(t,x,y) — E x [t > t d , q(t - t d , X Td , y)), x,y G D 

is continuous and the transition density of X D . For any bounded open set D C M. d , we will use Gd 
to denote the Green function of X D , i.e., 

/•OD 

G D (x,y) := q D {t,x,y)dt, x,y G D. 
Jo 

Note that Gd{x, y) continuous in (D x D)\{(x,x) : x G D}. We will frequently use the well-known 
fact that Go{-,y) is harmonic in D \ {y}, and regular harmonic in D \ B(y, e) for every e > 0. 
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Using the Levy system for X, we know that for every bounded open subset D, every / > and 
all x £ D, 

E x [f(X TD ); X TD _ £ X TD ] = [_ [ G D (x,z)J(z - y)dzf(y)dy. (4.16) 

Jd c Jd 

For notational convenience, we define 

K D (x,y) := f G D (x,z)J(z-y)dz, (x,y)eDxD c . (4.17) 
Jd 

Thus (|4.16p can be simply written as 

E x [f(X TD ); X TD _ t X TD ) = [_ K D (x,y)f(y)dy, (4.18) 

Jd c 

revealing K£>(x,y) as a density of the exit distribution of X from D. The function K£>(x,y) is 
called the Poisson kernel of X. Using the continuity of Go and J, one can easily check that is 
continuous on D x D C . 

The following proposition is an improvement of Lemma 14.31 The idea of the proof comes from 

EH. 



Proposition 4.9 For all r > and all xq E R rf , 

E x [t bm \ < 2V(2r)V(r-\x-x \), x E B(x Q ,r) . 
In particular, for any R > 0, r E (0, R) and xq E M. d , 



^x[t BM ] < C 7 (0(r- 2 )^((r-|x-x o |)- 2 ))- 1/2 

r a/2 ( r _ | a ._ a; |)a/2 

~ (£( r -2))l/2 (^( r _ Ix-Xol)- 2 )) 1 / 2 ' 

where Cj = Cj{R) is the constant form Proposition Iff. 131 



x E B(x Q ,r) , 



Proof. Without loss of generality, we may assume that x = 0. For x ^ 0, put Z t = 4±f . Then 
Zt is a Levy process on R with 

E(e iflZt ) = E(e i9 ^ Xt ) = e^H^ = e"*^ 2 ) 9 G R. 

Thus Z t is of the type of one-dimensional subordinate Brownian motion studied in Section 13.31 It 
is easy to see that, if X t E 5(0, r), then \Z t \ < r, hence 

E zb"B(0,r)] < ^\x\[t], 

where f = inf{i > : \Zt\ > r}. Now the desired conclusion follows easily from Proposition 13.131 
(more precisely, from (|3.28p ). □ 

As a consequence of Lemma 14.21 Proposition 14.91 and (|4.17p . we get the following result. 
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Proposition 4.10 There exist C15, C±q > such that for every r G (0, 1) and xq G R , 

K B{xo>r) (x,y) < C 15 j(\y-x \-r){<p(r- 2 )4>((r-\x-x \)- 2 ))- 1/2 (4.19) 

r a/2 ( r _ | x _ X() |)«/2 



j(|y - sol - r) 



(£( r -2))l/2 (^(( r _ (a- _ 3, |)-2))l/2 ' 



/or aZZ G B(xo,r) x B(xq, r) and 



K B(xo , r) (x ,y) > C M ^'y ?P xjXIy-xol)^ (4-20) 



/or a// y £ ^(^o, r) . 

Proof. Without loss of generality, we assume xq = 0. For z G 5(0, r) and r < |y| < 2 

|y| - r < \y\ ~ \A <\z-y\< \z\ + \y\ < r + |y| < 2|y|, 
and for z G 5(0, r) and y G 5(0, 2) c , 

\y\-r < \y\ - \z\ <\z-y\< \z\ + \y\ < r + |y| < |y| + 1. 
Thus by the monotonicity of j, f|3. 15|) and (|3.16p . there exists a constant c > such that 



cj(|y|) < j(\z-y\) < j(\y\-r), (z,y)eB(Q,r) x 5(0, r) . 

Applying the above inequality, Lemma 14.21 and Proposition 14.91 to (|4.17|) . we have proved the 
proposition. □ 

Proposition 4.11 For every a G (0, 1), r G (0, 1/4), xq G R d and x\,X2 G B(xq, ar), 



K B(x ,r)(xi,y) < C li K B ^ XOi ^(x 2 ,y), y£B(x ,r) 
where C\\ = C\±{a) is the constant from Corollary \4-8\ 



Proof. Let a G (0, 1), r G (0, 1/4) and x G R d be fixed. For every Borel set A C 5(xo,r) , the 
function x 1— > V X (X TB ^ r) G A) is harmonic in B(xQ,r). By Corollary 14.81 and (|4.18|) . we have for 
all x\,X2 G B(xo,ar), 

KB(x ,r)(xi,y)dy = F X1 (X TB G A) 

J A 

< C U F X2 {X r) e A) = / K BM (x 2 ,y)dy. 

J A 

c 

This implies that K B ^ xo r ^(xi,y) < CuK B f X0tr \(x 2 ,y) for a.e. y G B(xo,r) , and hence by the 

C 

continuity of K B ^ XQ ^(x, •) for every y G B(x , r) . □ 
The next inequalities will be used several times in the remainder of this paper. 
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Lemma 4.12 There exists C > such that 



s a/2 a/2 

Tjo < C TT? > < s < r < 4, (4.21) 

s l-o/2 r l-a/2 

Tjk < C -r, < s < r < 4, (4.22) 

(£(s- 2 )) 1/2 (£(r~ 2 )) 1/2 

, a ~ a / 2 (^(s- 2 )) 1/2 < Cr 1 -"/ 2 (^(r- 2 )) 1/2 , < s < r < 4, (4.23) 

(£(s- 2 )) 1/2 (^(r- 2 )) 1/2 , s 

(^( S -2))i/2 ^(r- 2 )) 1/2 t 

^*< C 4£>, 0<r<4, (4.26) 

^ S <C^5, 0<r<4, (4.27) 



o s 



and 



r s*- 1 „ r 



a 4 s") ds - c n^-y 0<r£i - (4 - 28) 

Proof. The first three inequalities follow easily from |3l Theorem 1.5.3], while the last five from 
the 0- version of (3] 1.5.11]. □ 

Proposition 4.13 For every a G (0, 1), there exists Cyj = Cn{a) > such that for every r £ (0, 1) 
and Jo € R d , 

r a/2-d ^((| y _ XQ | _ r )-2))l/2 



(^(r-2))l/2 (\y - XQ \ - r )a/2 

Vx 6 B(xq, ar), y £ {r < \xq — y\ < 2r} . 



Proof. By Proposition 14.111 

K B (x ,r)(x,y) < ^ / K Bix0tr) (w,y)dw 

' J B(xQ,ar) 

for some constant c\ = c\{a) > 0. Thus from Proposition 14.91 (|4.19p and Remark 13. 101 we have that 



K B(x ,r)(.x,y) < ■% / / G B r X0 Aw,z)J(z-y)dzdw 

r JB(x ,r) JB(x ,r) 



r d Jb(x ,v) 



^z[T B (x ,r)]J{z ~ V)dz 



c 2 r"/ 2 r (r-|z-x |)"/2 _ 

" (£(r-2))i/2 y B(xo r) (*(( r _ | z _ xo|)-2))V2 ^ y > Z 



35 



for some constant C2 = 02(a) > 0. Now applying Theorem 13.41 we get 



C3r a/2-d j- ( r _ | z _ Xo |)«/2 £{\z-y\- 2 ) 



K BM (x,y) < (£(r _ 2))1/2 {i{{r _ lz _ Xolr2)) i/2 \ z _ y \d + c 

for some constant C3 = 03(a) > 0. Since r — \z — xq\ < \y — z\ < 3r < 3, from (|4.2ip we see that 



(r-\z-x \) a / 2 . {\y-z\) a/2 



(£((r- |z-x |)- 2 ))V2 - ^(£(| y _ z |-2))i/2 
for some constant C4 > 0. Thus we have 

C5rQ/2 - rf /■ w , |- 2)) i/2 



(£(r-2))l/2 7 B(yi|y _ xo| _ r)c |z-y|^/2 
< c 6 r"/ 2 - rf /- (^' 2 )) 1/2 ,, 



for some constants C5 = 05(a) > and cq = c§(a) > 0. Using (|4.24j) in the above equation, we 
conclude that 

*B(x ,r){X,y) - Wr - 2)) l/2 { \ y _ XQ \ _ r)a /2 

for some constant C7 = 07(a) > 0. □ 



Remark 4.14 Note that the right-hand side of the estimate can be replaced by 
4.3 Boundary Harnack principle 



V(r) 



r d V(\y-x \-r) ' 



In this subsection, we additionally assume that a £ (0, 2 A d) and in the case d < 2, we further 
assume (|3.2|) . 

The proof of the boundary Harnack principle is basically the proof given in [25J, which is adapted 
from [HU2]. The following result is a generalization of |42l Lemma 3.3]. 

Lemma 4.15 For every a £ (0, 1), there exists a positive constant C±g = C\g(a) > such that for 
any r £ (0, 1) and any open set D with D C -6(0, r) we have 

F X (X TD G B(0,r) c ) < C 19 r- a l{r- 2 ) I G D (x,y)dy, x £ D n B(0, ar). 

Jd 

Proof. We will use C™(R d ) to denote the space of infinitely differentiable functions with compact 
supports. Recall that L is the L 2 -generator of X in (|4.ip and that G(x,y) and Go(x,y) are the 
Green functions of X in M. d and D respectively. We have LG(x,y) = —S x (y) in the weak sense. 
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Since Go{x^y) = G(x,y) — K x [G(X TD ,y)], we have, by the symmetry of L, for any x G D and any 
nonnegative (f> G C£°( 



G D {x,y)~L(f)(y)dy = / G D (x,y)~L(f)(y)dy 
l • J«. d 

G(x,y)L<Ky)dy - [ E x [G(X TD ,y)]L(j)(y)dy 

JR d 

G(x,y)L<l>(v)dy- f [ G{z,y)^{y)dy¥ x {X TD G dz) 
Jd c Jti d 

- ( f>(x) + E x [ ( p(X TD )}. 



= -<P(x)+ / cf>(z)¥ x (X TD £dz) 
In particular, if (f>(x) = for x 6 D, we have 



E X [4>(X TD )] = / G D (x,y)Lct>(y)dy. 

JD 



(4.29) 



For fixed a G (0, 1), take a sequence of radial functions 4>m in C c °°(R d ) such that < <j) m < 1, 



0, |y| < a 

1, 1 < |y| < m + 1 
0, |y|>m + 2, 



and that Yi,j 



dyidyj 



is uniformly bounded. Define 4> m ,r(y) = 4>m(-) so that < <j)m,r < 1> 



<t>m,r{y) = < 



'0, 1 2/| < ar 

1, t < \y\ < r(m + 1) and sup 
[0, |y|>r(m + 2), 



dyidyj 



■(v) 



< ci r 



We claim that there exists a constant ci = ci(a) > such that for all r G (0, 1), 

sup sup \L(j) m>r (y)\ < ar~ a £(r~ 2 ). 

m>l i 



(4.30) 



(4.31) 



In fact, by Proposition 13.41 we have 

[4>m,r{x + y) ~ 4>m,r(x) - {V (t>m,r{ X ) ' V) l B(Q,r){y))J {y)dy 



< 



{\v\<r} 



.{X + y)~ <prn,r( X ) ~ {^4>m,r(x) ■ y) l B (Q, r ) {y))J{y)dy 



+2 / J{y)dy 

J{r<\y\} 

<% [ \y\ 2 J(y)dy + 2 [ J(y)dy 

r J{\y\<r} J{r<\y\} 



<-%[ 



1 



{\y\<r} \V\ 

c,-2\ 



d+a-2 



l )dy + C 3 / 
Jh 



{r<\y\} 



\y\ d+ ' 



A\y\~ 2 )dy 



r 2 J s a 



—ds + C4 



l(s- 2 ) 



,1+Q 



ds. 
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Applying (|4.26p - (|4,27p to the above equation, we get 



(x + y) -4> m ,r{x) - (V<p mtr (x) ■ y)l B (0,r) {y))J{y)dy 



< c 5 r- a £(r- 2 ) 



for some constant C5 = 05(0*, a, •£) > 0. So the claim follows. Let A(x,a,b) := {y G M. d : a < 
\y — x\ < b}. When D C .6(0, r) for some r G (0, 1), we get, by combining (|4.29p and (|4.3ip . that 
for any x G D n B(0, ar), 



\(X TD GB(0,r) 



lim F X (X TD G A(0,r,(m + l)r)) 



m— >oo 
-a 



< C7r"^(r- 2 ) / G D (x,y)dy. 



D 



□ 



Lemma 4.16 There exists C20 > snc/i £/iat /or any open set D u>ii/i -B(j4, kt) d D <Z 13(0, r) for 
some r G (0, 1) and k G (0, 1), we have that for every x G D \ B(A, nr), 



G D (x,y)dy 
D 

nr\—2 



Proof. Fix a point x G D \ S(il,Kr) and let 5 := B(A, ^f). Note that, by the harmonicity of 
Gd(x, ■ ) in D \ {x} with respect to X, we have 

G D (x,A) > f K B (A,y)G D (x,y)dy > f K B (A,y)G D {x,y)dy. 

J DnB c JDr\B(A,2zL)c 

Since < \y — A\ <2r for y G B(A, 2ip) c p 2) anc | j j s a decreasing function, it follows from 
P~20j) in Proposition g3D] and Theorem E31 that 

> d ^ / G D (x,y)J(y-A)dy 



for some positive constants ci and C2. On the other hand, applying Theorem 14.71 we get 
/ G D (x,y)dy<c 3 G D (x, A)dy < c A r d K d G D (x, A), 
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for some positive constants C3 and C4. Combining these two estimates we get that 

?D(x,y)dy < c 5 (rV + r^ ^J^ G D (x,A) (4.32) 

for some constant C5 > 0. 

Let n = D\B(A,f). Note that for any z G B(A,f) and y G fi, 2~ 1 |y-z| < |y-A| < 2\y-z\. 
Thus we get from <@TQ and (f3TT5|) that for z G -B(A, 



c^lfafoA) < K u (x,z) < c 6 K n (x,A) 



(4.33) 



for some cq > 1. Using the harmonicity of Gd{-,A) in D \ {A} with respect to X, we can split 
Gd(-,A) into two parts: 



G D (x,A) =E X [G D (X Tn ,A)j 



+E, 



G D (X rn ,A): X Tn G{^<|y-^|<y} 



:= /1+/2. 

Since Gd{v,A) < G(y,A), by using (|4.33|) and Theorem 13.21 we have 



h < c 6 K n (x,A) / G D (y,A)dy 

< Cl K a(x ,A)f w - JF -L—^ ¥) 



dy, 



for some constant c-j > 0. Since \y — A\ < 4r < 4, by (|4,2ip . 

\y - A\ a / 2 Ur) a / 2 
< c 8 



" s ^((4r)-2) 
for some constant eg > 0. Thus 

ii < c 7 c 8 ^(x,A) / 
Je 



(4r 



, a/2 



|y-A|^/2£((4r)-2) 



1 



^((4r)- 2 7 

for some constant eg > 0. Now using (|4.33p again, we get 

1 



h < c 10 K a ' z - a r 



a/2—da~d_ 



c W K a / 2 - d r a - d 



1 



£((4r)- 2 ) 



2\ u * I Ti 



(4.34) 
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for some constant c\q > 0. On the other hand, again by Theorem 13.21 and (|4.34[) 

h = I G D (y,A)F x (X Tn £dy) 
J{^<\y-A\<^} 

5 c " 4a»-^?> <(I»-V») p ' (jr '» 6 *» 



for some constants en > and c\i > 0. 
Therefore 



1 _ „, . KT, 



for some constant C13 > 0. Combining the above with (|4.32p . we get 

G D {x,y)dy 

D 

for some constant C14 > 0. It follows immediately that 

/ G D (x,y)dy 
Jd 

< r« (l + ||^|) P, (X TD , W G B(A, «r) 

□ 

Combining Lemmas 14. 15ll4TT6l and using the translation invariant property, we have the following 

Lemma 4.17 There exists C21 > such that for any open set D with B(A,nr) C D C B(Q,r) 
for some r £ (0, 1) and k G (0, 1), we have that for every x £ D n -B(Q, §), 

P,(Jr TD e B(Q,r)°) 

Let a, b) := {y £ R d : a < \y - x\ < b}. 



Lemma 4.18 Let D be an open set and r 6 (0, 1/2). For every Q G M rf and any positive function 

iRr, -ir) suc/i t/iai /or any xeflfl -B(Q, §? 



n vanishing on D c n if r )> ^ere is a a £ (-jpf", -ir) suc/i t/iai /or any ifDn |r), 



E, 



«(^ DnB(Q , CT) );X TDnS{Q CT) G 5(Q,a) c ] < C 22 jf ^ ^ J(y - Q)u{y)dy (4.35) 



/or some constant C22 > independent of Q and u. 
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Proof. Without loss of generality, we may assume that Q = 0. Note that by (|4.25p 

£((\y\-a)- 2 )^ 2 (\y\-a)-^ 2 u(y)dyda 



f-r JA(Q,a,2r) 



A(0,fr,2r) J ±gr 



i{{\y\-a)- 2 ) ll2 (\y\-o)- al2 dou(y)dy 



< 



A(0,^r,2r) Wo 



< Cl 



A(0,^,2r) 



2\ X /2 / ln x 



\y\- — 



lOr V 



for some positive constant c\. Using (|4.22p and (|4.23p . we get that there are constants ci > and 
C3 > such that 



\y\ 



10r 
~6~ 



1/2 



10r y-«/2 
|y| - -g- J 



< c 3 / l{\y\~ 2 ) l l 2 \y\ l - a l 2 u(y)dy 

£((2r) ^J 1 / 2 yA(o,ifc,2r) 
Thus, by taking C4 > 6C1C3, we can conclude that there is a a 6 (^§ r , ir r ) such that 



/ i((\y\-aT 2 ) 1/2 (\y\-vr a/2 <y)dy 

JA(0,a,2r) 

l{\y\~ 2 )u(y)dy. 



.-a/2 



< C4 



£((2f)-2)l/2 y j4(0i lor^ 



(4.36) 



Let 1 £ Dfl 5(0, |r). Note that, since X satisfies the hypothesis H in [43J, by Theorem 1 in 
we have 



u (^DnB(o, CT )); x T Dns(0iCT) g 5(0, a; 

u ( 1 ^nB(0,,));^Dn B( o ! a) G 5(0, cr) 6 , T DaB(0i(j) = r B(0)O -) 
U ( X t B (0,<t))i X -r B (0,a) G - B (0 ) <T ) C > r flnB(0, ff ) = T B(0,a) 



^ X r B(0 ^X TB(0 ^ €5(0, a) 



B(0,cr) c 



K B(o,a)( x ,y)u(y)dy. 
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Since a < 2r < 1, from (|4.19p in Proposition 14. 101 Proposition 14. 131 we have 



u (XT DnB ( () X TDnB G B(0,a) c < / K B{0 ^ ) {x,y)u{y)dy 

'B(0, CT ) C 



r W 2 - rf {i{{\y\-v)- 2 )) 1/2 ( w 

■/B(0,2r)<= (^(O- 2 )) 1/2 (£((<?- Pi) 2 )) 1/2 

for some constant C5 > 0. 

When y E ^4(0, 2r, 4) we have j%\y\ < |y| — <r, while when |y| > 4 we have |y| — a > \y\ — 1. 
Since a — |x| < <r < 2r, we have by (|4.2ip and the monotonicity of j, 

m - a) (Ipp - |s|)-*))V2 ^ C6J (,12 J ^(27p) ' ^^°' 2r ' 4 ) 



and 

W 2 (g - [xl) a / 2 

'(^(ct- 2 )) 1 ^ (l{(o- \x\)- 2 )) 1 / 2 ~ ^ JVIWI ^((2r) 

for some constant C6 > 0. Thus by applying (|3. 15 j) and (|3.16p . we get 



idyl - ttt^t^t^ ,'_,w^m/9 < c 6j(l?/l-i) „^_2v 4 



fJ Q / 2 (cr-ldW 2 



for some constant cj > 0. Therefore, 



JB(0,2r) c 



W 2 (a-|x|) Q / 2 
(%- 2 ))V2 (^((a-lxD-a))!^"^ 1 ' 



On the other hand, by (|4.2ip . (|4.36p and Theorem 13.41 we have that 

JA(0,«,2r)(£(a-*)f 2 (\v\-<r)*' 2 U[y)dy 

(10rY d W 2 f ffl(|y|-a)- 2 ))V 2 
" ItJ (iiopC, (|y|-.)"/ 2 U(y)dy 



< c 8 r~ d — — / £(|y|- 2 )n(y)dy 

W(2r)- 2 )) 1 / 2 (^((2r)- 2 )) 1 / 2 A(o,±te 2r) U ' ' K ' 



^ °» p«l\-2^ I , m- 2 )\y\- d ' a u(y)dy 

< c io ^/ w9 . / J(y)u{y)dy 



7((2r)" 
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for some positive constants eg, eg and cio- Hence, by combining the last two displays we arrive at 



U ( X r D n B(0 ,^ X r D nB^) G B ^^T\ < C H l(( 2 r)~2) / fl Q jg, J ivXv)^ 



'B(0,i§E)«= 

for some constant en > 0. □ 

Lemma 4.19 Lei I? 6e an open se£ and r G (0,1/2). Assume that B(A,kv) C D (1 B(Q,r) for 
k G (0, 1/2]. Suppose that u >0 is regular harmonic in D n B(Q,2r) with respect to X and u = 
in .D c n .B(Q, 2r). If w is a regular harmonic function with respect to X in D n r) suc/i i/iai 



i£B(Q,|) c Up c nB(g,r)), 



0, xGi^r'l), 



.^((2r)~ 2 ) ^ _ ^ n ^ D ^ 3 

/or some constant C23 > 0. 



> w(A) > C 23 « a „))J' n(x), xefln B(Q, -r) 



Proof. Without loss of generality we may assume Q = 0. Let x G Dn5(0, |r). The left hand side 
inequality in the conclusion of the lemma is clear from the fact that u dominates w on (Df]B(0, r)) c 
and both functions are regular harmonic in D(~)B(0, r). Thus we only need to prove the right hand 
side inequality. By Lemma [4. 181 there exists a G -g^-) such that (|4,35p holds. Since u is regular 
harmonic in D n 5(0, 2r) with respect to X and equal to zero on D c n -B(0, 2r), it follows that 

u(x) = E x [u(X TDnBi0 ^); X TDnB{0 a) G B(0,a) c \ < ci^py J b{o m)e J(vHy)dy (4.37) 
for some constant c\ > 0. On the other hand, by (|4.20p in Proposition 14.101 we have that 
'{A) = / K DnB{0r) (A,y)u(y)dy> K B(A ) (A,y)u{y)dy 



w{ 

(nr) ( 



> c 2 I J(A-y) 1 ±-L— U {y)dy 



for some constant C2 > 0. Note that \y — A\ < 2\y\ in A(0, ^-,4) and that \y — A\ < \y\ + 1 for 
|y| > 4. Hence by the monotonicity of j, (|3.15p and (|3.16p . 



B(0,^) c 



for some constant C3 > 0. Therefore, by (|4.37p 



□ 
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Definition 4.20 Let k G (0, 1/2]. VFe say £aai an open sei -D in M. d is K-fat if there exists R > 
suc/i that for each Q G <9Z? and r G (0, R), D n B(Q,r) contains a ball B(A r (Q), nr). The pair 
(R,k) is called the characteristics of the K-fat open set D. 

Note that all Lipschitz domain and all non-tangentially accessible domain (see |24j for the 
definition) are K-fat. The boundary of a K-fat open set can be highly nonrectifiable and, in general, 
no regularity of its boundary can be inferred. Bounded K-fat open set may be disconnected. 

Since £ is slowly varying at oo, we get the Carleson's estimate from Lemma 14.191 

Corollary 4.21 Suppose that D is a K-fat open set with the characteristics (R,k). There exists 
a constant C24 depending on the characteristics (R,k) such that ifr<RA^,Q^ dD, u > is 
regular harmonic in D n B{Q, 2r) with respect to X and u = in D c n B(Q, 2r), then 

u (MQ)) > C24 u(x) , Vx G D n B(Q, ^r) . 

The next theorem is a boundary Harnack principle for (possibly unbounded) K-fat open set and 
it is the main result of this subsection. 

Theorem 4.22 Suppose that D is a K-fat open set with the characteristics (R,k). There exists a 
constant C25 > 1 depending on the characteristics (R, k) such that if r < R A I and Q G dD, then 
for any nonnegative functions u,v in W 1 which are regular harmonic in D n B(Q,2r) with respect 
to X and vanish in D c n B(Q, 2r), we have 

r -i <MQ)) < u{x) u(A r (Q)) rr ^ nnR(n r \ 

Proof. Since I is slowly varying at 00 and locally bounded above and below by positive constants, 
there exists a constant c > such that for every r G (0, 1/4), 

( ^)~ 2 ) ^f)' 2 ) ^r 2 ) \ < c f438) 

V^M" 2 )' ^((4r)~ 2 )' ^((4r)- 2 )' ^((2r)- 2 ) J ~ { ' 

Fix r G (0, R A 4) throughout this proof. Without loss of generality we may assume that Q = 
and u(A r (0)) = v(A r (0)). For simplicity, we will write ^4 r (0) as A in the remainder of this proof. 
Define u\ and ui to be regular harmonic functions in D n -B (0, r) with respect to X such that 



ui{x) 

and 



u(x), ieA(0,r,|), 

0, x G B(0, ^) c U (D c n B(0, r)) 



0, xGA(0,r,f), 



n(x), x G -B(0, 4f ) c U (D c n B(0, r)). 
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and note that u = u\ + u%. HDD A(0, r, ^) = 0, then u\ = and the inequality (|4.42p below holds 
trivially. So we assume that D n A(0, r, ^) is not empty. Then by Lemma 14.191 



u(y) < c x k 



AW 



for some constant c\ > 0. For x G D D -8(0, |), we have 



3r 



3r , 



T I)nB(0,r) , 



X r D n B(0 , r) €DnA(0,r,- 



< [ sup u{y) ] P x . ( X TDnB{0 r) G £> n A(0, 



v OnA(0,r,|) 



3r 



< ( sup n(y)]P x .(X TcnB(Or) Gi?(0,r 



v DnA(o,r,|:) 



< Cl k' 



u(A)F x [X TDnB(o .eB(0,r) 



Now using Lemma liTfl fwith D replaced by Dt~)B(0, r)) and (|4.38p . we have that for x G DnS(0, §), 



Ml(x) 
< C 2 K 



(4.39) 



rf - § « ^((^)- 2 ) ^- 2 ) l , wrn uU) x 

£((2r)- 2 ) £((4r)- 2 ) V ^((4r)~ 2 ) y 1 ; 



X ^ l^(DnB(0,r))\B(A,f) G B {A, 

< c 3 ^)P x ^ T(DnB(0 r)AB(A ^ G B(A -) 



T 



(4.40) 



for some positive constants c 2 and C3 = C3(k). Since r < 1/4, Theorem 14.71 implies that 



u{y) > c A u(A), y G B{A,—) 



for some constant C4 > 0. Therefore for x G D H 5(0, |) 

u(s) = E x . [«(X r(J3nfl(0jr))Xfl(A iJf) )] > c 4 u(A)F x (x T(DnB{0r 



nr . 



(4.41) 



Using (|4.40p . the analogue of ()4.4ip for v, and the assumption that u(A) = v(A), we get that for 
x G flnB(0, |), 



«i(x) < ca^P^X^^^^ 6fl(A ., 
for some constant C5 = cs(k) > 0. For x £ D (1 B(0, r), we have 



< c 5 v(x) 



(4.42) 



u 2 (x) 



K D nB(o,r)(x,z)u(z)dz 



B(0,^-Y JDnB(0,r) 



GDnB(o,r) (x, y)J(y- z)dy u{z) dz. 
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Let 



s(x) := G DnB ^ r) (x,y)dy. 

JDnB(0,r) 



Note that for every y G 5(0, r) and z G B(0, ^) c , 

1. , , , , , , , . 

gp| < p| — r < p| — 1 2/ 1 < \y — z \ < |y| + I- 2 ! < r+ p| < 2|z| 

and that for every y € B(0,r) and z G B(0, 12) c , 

1 2r| — 1 < \y — z\ < \z\ + 1. 

So by the monotonicity of j, for every y G -B(0, r) and z G ^4(0, , 12), 

j(12M) < j(2\z\) < J(y -z) < j f i|H < j f M 

and for every y G 5(0, r) and every z G B(0, 12) c , 

i(N-i) < < i(N + i)- 

Using (|3.15p and (|3. 16[) . we have that, for every y G -6(0, r) and z G -B(0, ^) c , 

c^VCkl) < J(v-z) < cei(N) 

for some constant C6 > 0. Thus we have 

for some constant cj > 1. Applying (|4.43j> to u, and t> and Lemma 14.191 to v and V2, we obtain for 
x G On5(0, |), 

/ x / ^ , u 2 (i) , , n £((Kr)~ 2 )u(A) . . 

s(A) v 2 (A) £{{2r) z ) v(A) 

= c8K ~ a jmSj v2ix ^ (4 - 44) 

for some constant c$ > 0. Combining (|4.42p and (|4.44p and applying (|4.38p . we have 

u(x) < c 9 v(x), x G D n -B(0, -), 
for some constant eg = cg(/c) > 0. □ 

Acknowledgment: We thank Qiang Zeng for his comments on the first version of this paper. We 
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